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^ ; ABSTRACT 

We present an analysis of 7 clusters observed by XMM-Newton as part of our survey of 17 most X-ray 
^ luminous clusters of galaxies at 2; ~ 0.2 selected for a comprehensive and unbiased study of the mass 
Q\ ■ distribution in massive clusters. Using the public software FTOOLS and XMMSAS we have set up an 
O , automated pipeline to reduce the EPIC MOS & pn spectro-imaging data, optimized for extended sources 
I analysis. We also developped a code to perform intensive spectral and imaging analysis particularly focussing 
on proper background estimate and removal. XMM-Newton deep spectro-imaging of these clusters allowed 
O ■ us to fit a standard /?-model to their gas emission profiles as well as a standard MEKAL emission model to 
their extracted spectra, and test their inferred characteristics against already calibrated relations. 
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INTRODUCTION 

Being the most massive gravitationally bound objects in the Universe, clusters of galaxies are prime 
targets for studies of structure formation and evolution. 
5^ , Of particular interest are the most X-ray luminous clusters. Indeed, it has been convincingly demonstrated 
that high X-ray luminosity is a reliable indicator of cluster mass {e.g. [Allen, 1998]). Therefore, X-ray 
selected clusters samples {e.g. [Ebeling et al., 1996]) are the best choice for undertaking representative 
studies of the mass distribution in massive. X-ray luminous clusters. In particular, as they are easily 
detected out to very large redshift, we can hope to probe cluster evolution back from early epochs. However, 
X-ray luminosity is also sensitive to the presence of cooling flows, mergers and non-thermal effects, thus 
the Luminosity-Mass and Temperature-Mass relation needs to be clearly understood before any strong 
ascertion is made on cluster evolution ([Smith et al., 2003]). Proper mass estimate is likely to arise only by 
combining multiwavelength observations of clusters of galaxies: namely strong and weak lensing observations 
([Smith et al., 2001]), velocity distribution and dynamics of the cluster galaxies (in a similar way as in 
[Czoske et al., 2001]) and through accurate, spatially resolved X-ray temperature measurements as we are 
concentrating on here. 

Following this premise, we are conducting a survey (Table 1) of the most X-ray luminous clusters 
(Ljffo i_2 4 kev] — ^ ^^^^ ^''^9 ' ^^^) ™ ^ narrow redshift slice at 2; ~ 0.2, selected from the XBACS (X- 
ray Brightest Abell-type Clusters Sample; [Ebeling et al., 1996]). As XBACS is restricted to Abell clusters 
([Abell, Corwin &; Olowin, 1989]), it is X-ray flux limited but not truly X-ray selected. However, a com- 
parison with the X-ray selected ROSAT BCS (Brightest Cluster Sample, [Ebeling et al., 2000]) shows that 
~ 75% of the BCS clusters in the redshift and X-ray luminosity range of our sample are in fact Abell clusters. 
Hence, our XBACS subsample is, in all practical aspects, indistinguishable from an X-ray selected sample. 
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Table 1. Our cluster list sorted by right ascension. Redshift, Hydrogen column density, flux and luminosity 
(columns 3 to 6 resp.) are quoted from the XBACS updated database. Out of these 17 clusters, 8 have 
HST/WFPC2 observations performed (as part of cycle 8 [ID:8249, PhKneib]); 3 only have a shallow HST 
observation through a snapshot survey (8301 &i 8719 programmes, PhEdge); 2 have not yet been observed 
with HST; the remaining 4 having archival data. Moreover, 11 have been observed in 1999/2000 with the 
wide field CFH12K camera ([Czoske et al., 2003]). The last column shows the effective exposure duration 
(after flares removal) for clusters already observed by XMM-Newton (resp. MOS and pn instruments), or 
quote the priority level within our proposal for those not yet observed. 
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THE ANALYSIS PIPELINE 

As to the version 5.2 of the official XMM-Newton data analysis software {XMMSAS)^ the extended 
sources could not be handled properly, regarding the background subtraction and vignetting corrections. 
The development of methods, where weights are applied to events according to their energies and positions 
as well as calibration data so as to correct the event list as if the instrument response were uniformely 
identical to the on-axis response, started as soon as 2001 ([Arnaud et al., 2001, Majerowicz et al., 2002, 
Marty et al., 2002, Marty et al., 2003]). We concentrated on our side to embedding these algorithms into 
an environment capable of pipeline processing a series of data from clusters of galaxies observations. Indeed, 
the manual operation of XMMSAS (using either external correction routines or integrated algorithms from 
versions 5.3 and later) is extremely time-consuming: several days for each dataset. 

Thanks to our pipeline, we reduced the required time down to less than 10 h per dataset, adding a working 
day for eye-inspection of the full sample results and another for compilation and final manual calculations. 

This pipeline integrates the latest calibration data from the EPIC consortium as well as specific fitting 
routines written in IDL and XSPEC environments. It is fully configurable to allow eye-inspection and 
partial or total re-runs if needed. It is extensively detailed by [Marty et al., 2003], but useful information 
about "blank" and "dark" fields (for background estimation and subtraction) are also respectively given by 
[Lumb et al., 2002, Marty et al., 2002]. The general processing scheme is the following: 

• re-process each cluster raw data to account for latest calibration knowledge; 

• re-project "blank" and/or "dark" data onto telescope aspect solution used for each cluster; 

• extract valid events within good-time- intervals (to reject noise and proton fiares); 

• detect and mask out spurious components (field sources and remaining bad pixels); 

• weight the events according to calibration data; 

• find the best cluster centroid; 

• process radial brightness profiles and fit a standard /3-model convolved by the instrument PSF; 

• process full maps in physical fiux units and adaptively smooth according to PSF description; 

• process spectra and call XSPEC for fitting (global for the mean temperature, annuli for the profile); 

• compute hardness maps and derive another temperature profile. 
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THE PIPELINE RESULTS 
Surface brightness profiles 

The core radius Rc, power index f3, normalization /o and background constant BKG are coming from 
the standard /3-model fitting. Note that the surface brightness profiles have been extracted between 0.3 and 
6.3 keV, logarithmically rebinned to enhance their signal-to-noise ratios and that the inner 40 arcsec of 
each cluster profile have been ignored for the fit so as to avoid fitting any emission excess, associated for 
example to a cooling flow. The background subtraction concerned only the instrumental component so that 
the BKG constant is a measure of the sky X-ray emission background. The /5-models are convolved by the 
instrument PSF before fitting. 

The Ri„ parameter is a measure of the radius inside which the cluster signal-to-noise ratio is greater 
than 1. The Fx {data) and Fx (model) parameters then are the integrated fluxes within this detection limit 
radius, from the data and from the best fit /?-model; in each case, these fluxes are cleaned from the sky 
background component; the difference between them thus shows the contribution of any excess emission. 

Surface brightness maps 

Surface brightness maps that have been produced in the same energy band during the pipeline analysis 
are not presented herein, because of room restrictions and because they are not necessary for the following 
calculations. However, they may be found on the internet^ if so desired. 

Spectral analysis 

Spectra are extracted in the 0.2 to 12.0 keV band and rebinned according the Poisson statistics require- 
ments (at least 25 net counts per bin). They are background subtracted from both instrumental components 
(estimated using the "closed" data) and sky components (from a spectrum extracted in a peripheral annu- 
lus of each dataset, between the 600 and 800 arcsec radii) and then fitted against the TBabs x MEKAL 
photo-absorption and bremsstrahlung models. 

The redshifts z and hydrogen column densities Nh have been respectively fixed to the same optical redshift 
references as used in the original XBACS database, and to the weighted means of the Nh measurements 
from the [Dickey &l Lockman, 1990] database included within 1 deg of the cluster center. The resulting 
temperatures kT, abundance relative to solar composition /x and normalization Sq are then reported. These 
measurements come from spectra extracted inside a circle of 400 arcsec radius, and hence integrate any 
central excess emission; that is why temperatures may appear a bit low. 

Direct temperature profiles 

Spectra have also been extracted within 40 arcsec-wide concentric annuli following the same algorithm 
than for global spectra. Fitting them using the same XSPEC settings hence led to a direct measure of the 
temperature profile of each cluster. However, no PSF correction have been applied from one annulus to 
another thus probably leading to a general underestimation in the fitted temperatures, specially toward the 
center ([Markcvitch, 2002]). 

A last spectrum has been extracted and fitted, in the same region than that used for the /3-model fitting, 
i.e. an annulus starting at 40 arcsec and extending upto the arbitrary limit of 400 arcsec. This allows 
to compare between the global temperature with the one found outside the core region, which may host a 
cooling flow. 

Intrinsic parameters calculation 

The luminosity distance and 1 arcsec equivalent length at the cluster position are deduced from 
the cluster redshift z and the default cosmology (ilo = 1, f^A = and Hq = 50 /150 km/ s/Mpc). 

The virial radius Ry has been identified with the distance r2oo from the cluster center where the mean 
enclosed overdensity equals 200, as calibrated upon the measured X-ray temperature by [Evrard et al., 1996]. 

The X-ray luminosity Lx was integrated between our XMM-Newton energy bandwith (0.3 — 6.3 keV) 
within the virial radius defined above. We applied a numerical bolometric correction to deduce the absolute 
luminosity Lx{bol). 

'^Only compiled on ftp://www-statiorL.ias.u-psud..fr/pub/epic/HiLx/2002Marty_saiitiago_poster. jpg at the moment, but 
individual maps will soon be available in the same ftp directory. 
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The galaxies velocities dispersions a^ij^ and a^ig have been reprinted from two weak lensing analy- 
sis, respectively from NOT/ALFOSC and UH8k detectors ([Dahle et al., 1990]), and CFHT12k instrument 
([Bardeau & Kneib, 2003]). The corresponding /3 parameters were computed according to the relationship 
with the measured X-ray temperature in the isothermal hydrostatic equilibrium case ([Cavaliere & Fusco-Femiano, 1976]) 

The electronic density ne(0) at the center of the cluster has been estimated using the standard isothermal 
/3-model equations ([Cavaliere & Fusco-Femiano, 1976]) and the previously measured X-ray temperature and 
luminosity. By integration within the virial radius defined above, we also estimated the cluster gas mass. 
Finally, the total virial mass has been deduced from two different calibrated relations upon the measured 
X-ray temperature {Mb and My rcsp. from [Cavaliere & Fusco-Femiano, 1976, Arnaud & Evrard, 1999]) 
as well as from the virial theorem relation upon the galaxies velocities dispersion (My from the a^ig defined 
above, except for A773 for which we used a^ij^). The gas fraction fg simply is the ratio of the gas mass to 
the total mass (we chose Mb as the isothermal reference). 

DISCUSSION 

All analysed clusters show an excess in central brightness that seems correlated with the presence of 
a very bright central galaxy in optical images. They hence are all candidates for cooling flow hosts, but 
A267 and A383 temperature profiles do not exhibit significant cooling cores. In addition, all clusters but 
A1835 tend to follow the universal temperature profile calibrated by [Loken et al., 2002] (Figure 3), radii 
being rescaled either to core radius or virial radius natural units, with a declining slope toward the edges. 
They are also very close to the [DeCrandi & Molendi, 2002] polytropic interpretation of their Beppo-SAX 
observed clusters data, which show slopes more horizontal at high radii for those hosting a cooling flow than 
for non-cooling flow specimens, although our sample does not present a steeply cooling core in average. 

If the problem of the cooling core may be bypassed by using only outer regions for temperature and 
/?-model fitting, these declining profile still questions the hypothesis of isothermality followed throughout 
this first analysis. Nevertheless, we expect from the general good agreement of the different mass (Figure 2) 
and P parameters estimations that our optical {CFHT12k) and X-rays {XMM-Newton) measurements are 
still both quite consistent with the isothermal model. But we note also that all calculations rely on a virial 
radius that has been estimated using a relation based on the isothermal hypothesis; that illustrates the need 
for another independant estimator of this radius. Even the clue given by the temperature profiles (Figure 3), 
that seem to indicate that the virial radius is about 10 times the core radius, is not so helpful if models 
fitted to either optical or X-ray data still rely upon isothermality. 

We plotted a temperature-luminosity logarithmic diagram (Figure 1) for our subsample and compared 
it to the model proposed by [Arnaud &: Evrard, 1999]. We observe that five clusters are indeed well co- 
aligned but following a shallower slope. A209 seems to depart from this alignment, as well as A1835 in a 
very significant manner. We also compared a temperature-velocity dispersion diagram with the model from 
[Girardi et al., 1996] and found it to be very close to the lower 90% confidence level boundary, with again 
the exception of A1835. 

Finally, the temperature-mass diagram is remarkably consistent with a mean power index of 1.576, while 
the gas mass fraction seems to average about 10% except again for A1835, as well as A209 and A773. 

CONCLUSIONS & PERSPECTIVES 

While early XMM-Newton observations ([Arnaud et al., 2001, Majerowicz et al., 2002]) did not confirm 
declining temperature profiles as seen by late ASCA observations ([Markevitch et al., 1998]) and Beppo-SAX 
([DeCrandi & Molendi, 2002, Loken et al., 2002]), we find in the present analysis of 7 XBACS clusters that 
6 of them do follow this kind of polytropic model. The seventh, A1835, systematically departs from the 
average in all diagrams (temperature-mass, temperature-luminosity), which confirms the rather complicated 
nature of its internal dynamics ([Majerowicz et al., 2002]). 

However, we still are working toward the comparison of these emission weighted temperature profiles, that 
have been obtained without correction for PSF and may thus have been underestimated ([Markevitch, 2002]), 
with profiles built on the basis of hardness ratio maps, that may be adaptively handled to account for the 
PSF and binned between concentric isophotes rather than concentric annulii to adapt to non-spherical 
morphologies. Also, more work is intended so as to derive accurate error bars on each diagram, and attempt 
to estimate the virial radii from independent optical datasets. 
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Fig. 1. Left: temperature-luminosity diagram as compared to the 90% confidence level range from 
[Arnaud & Evrard, 1999] model. Right: temperature-velocity dispersion diagram as compared to the 90% confi- 
dence level range from [Girardi et al., 1996] model (diamonds show values from [Bardeau & Kneib, 2003], trian- 
gles from [Dahle et al., 1990]). 
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Fig. 2. Left: temperature-mass diagram, for different mass estimators, showing a mean power index of 1.576. 
Right: gas mass fraction, for the same estimators, as a function of the X-ray temperature. 
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Fig. 3. Temperature profiles with radius rescaled either to core radius or virial radius units. The dashed line 
stands for the average slope in the central part of the polytropic simulation of [Loken et al., 2002]. 
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Finally, this work has been made possible through the use of a pipeline analysis environment, specifically 
aimed at dealing with extended sources data from XMM-Newton instruments. It has confirmed the impor- 
tance of maintenance and upgrades, to keep up-to-date with the instruments calibration database, as well 
as the need for such automatic tools in the purpose of building catalogs from the latest X-ray telescopes, so 
that massive cluster investigation may be possible in the future: access to large statistical samples, cross- 
compare mass and morphological measurements with optical data, correlate X-ray and Sunyaev-Zel'Dovich 
data. . . The ESA's Planck satellite is promising thousands of millimetric clusters detections by the end of 
the decade, which could lead to narrower constraints on cosmology, provided an extensive X-rays database, 
like BAX ([Sadat et al., 2002]), also exists. 
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